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Fundamental problems of periodicity and transient process to periodicity of chaotic tra- 
jectories in computer realization with finite computation precision is investigated by 
taking single and coupled Logistic maps as examples. Empirical power law relations of 
the period and transient iterations with the computation precisions and the sizes of cou- 
pled systems are obtained. For each computation we always find, by randomly choosing 
initial conditions, a single dominant periodic trajectory which is realized with major por- 
tion of probability. These understandings are useful for possible applications of chaos, 
e.g., chaotic cryptography in secure communication. 
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1. introduction 

Chaos study has attracted much attention in the last half of the 20th century [1- 
7]. However, there have been still some fundamental problems in the chaos field 
remaining not completely understood, including the computer realization of chaos. 
A computer has always finite precision, and trajectories of any autonomous or pe- 
riodically driven deterministic systems must be periodic. Thus, we are faced with 
a difficult circle: chaos, which generally cannot be analytically shown, has to be 
manifested by numerical computations while any computer simulation can never 
produce a true chaotic orbit. 

In early time of chaos study, the commonly accepted way out of this difficulty 
is that the average period of computer realization of chaotic trajectories(CRCTs) 
is about the order of the digital space determined by the computer precision. For 
high precision (like double precision 2 -52 ) this period is very large (may be even 
unreachable in numerical simulations), and chaos can be clearly seen in the time 
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region much smaller than this period. Moreover, it is anticipated that from different 
initial conditions chaotic trajectories of a given dynamics system may eventually 
enter different periodic orbits, and the number of periods may be huge since there are 
infinitely many unstable periodic orbits embedded in any given chaotic attractors. 

However, the above intuitive understandings are not correct. In [8-11], the au- 
thors showed that the average length of the period with computation precision e 
is scaled as e~~ with D2 being the correlation dimension of the given chaotic at- 
tractor, which is much shorter than the possible length allowed by the computer 
precision. This short periodicity seriously influences some applications of chaos. 
For instance, chaos encryption has become an attracting topic in the recent ten 
decades [12- 18]. The authors in [19] suggested to use chaotic trajectories of a sin- 
gle map to reach extremely high level of security, but it was soon found [20] that 
the security is very low because the generated chaotic trajectories by the computer 
realization are periodic with rather short average period. 

Therefore, the problem of periodicity of CRCTs has both theoretical signifi- 
cance and practical importance, deserving further investigation. In this paper we 
will study this problem by taking single and coupled map lattices as our examples. 
We are particularly interested in the characteristic features of CRCT affecting the 
properties of chaotic cryptograph. In Sec. 2 we focus on the periodicity and the re- 
lated transient of a single map for different precisions, and find that both the average 
periods of CRCTs and the transients are indeed much shorter than the size of the 
digital space of computers in all the precisions tested, and these observations well 
agree with the theoretical predictions of [8,11]. A surprising observation is that for 
different initial conditions the number of periodic orbits of CRCTs is very small, and 
in all precisions tested one can find a single dominant period for randomly choosing 
initial conditions which is realized with a major portion of probability. In Sec. 3 we 
directly compute the period length of single logistic and tent maps by double preci- 
sion realization. In Sec. 4 we investigate the same problem for coupled systems, and 
find that the average period and transient iterations increases exponentially as the 
system size increases. This observation is very useful in some practical situations 
(such as chaos communications) where long periods are required. And the existence 
of dominant period is observed as well for coupled systems with different sizes. The 
last section gives brief discussion and conclusion where the significance of the results 
obtained in this paper for the applications of chaotic cryptography is emphasized. 

2. Periodicity and transient of a single chaotic map in finite 
precision computer realization 

We first consider the single logistic map as our example 

x n+1 = 4x„(l - x n ) (1) 

which is well known chaotic, and has aperiodic trajectory for an arbitrarily chosen 
typical initial xq in (0, 1) region. However, with finite-precision all trajectories must 
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be periodic, and the period of the finally realized motion depends on the initial 
condition. Our task in this section is to study how the average period and the average 
length of transient iterations towards periodic motions are related to the computer 
precision; and how many such periods can be found. All numerical investigations 
use Intel Pentium III 700MHz and Pentium IV 1.7GHz computers with Compaq 
Visual Fortran code. 

Let us follow the computer computation process with round-off truncations in 
finite precision. Suppose we take e = 10~ h (h = 1,2,...) computation precision, 
then for any given continuous x n value we measure only discrete value x n , with x n 
being an integer multiplier of e as 



In Eqs.(3) with the discrete operation the state x n — may be realized with certain 
finite probability [in fact this probability is zero for continuous variables, i.e., x n — 
cannot be realized unless xo takes certain nontypical values of zero measures]. 
Throughout the paper, we will exclude the x n = state solution from our counting. 
Map (3) is always periodic though map (1) is chaotic and aperiodic. As an example 
we consider h = 4, and start from x = 0.951636985290801. The iterations of (3) 
give 

x = 0.951636985290801 (x = 0.9516) -» zi = 0.18422976 (x 1= 0.1842) > 

x u = 0.453879 {x u = 0.4539) -> x 15 = 0.99149916 (x 15 = 0.9915) -> > (4) 

x 117 = 0.54599356 {x ur = 0.5460) -> i U8 = 0.991536 {x ns = 0.9915) 

For this initial condition the period of map (3) is T — 118 — 15 = 103, and the 
transient time for the trajectory to enter the periodic circle is r = 15. Of course, 
for different x 's we may find different periods and different lengths of transient 
iterations. 

With the above approach we can systematically investigate the dependence of 
t, T on the precision e = 10~ h . In Fig.l we plot r (circle) and T (square) vs h, 
where r and T are obtained as the averages with arbitrarily chosen different initial 
xq values, 




(2) 



With (2) the actual map becomes 



x n+1 = 4 x n (1- x n ) 
x n — _P(x n ) 



(3) 




(5) 
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where Tj and are the transient and period iterations computed by ith ini- 
tial condition, respectively, and i runs over all possible initial states for ft < 6 
(M = 10 h ) in the discrete variable space, and i is chosen randomly in (1, 10 h ) for 

/ M 

ft > 6. In Fig.l vertical bars indicate error estimates a T = J jj J2( T i ~ T ) 2 -> a T — 

I M 

\ ~M — T) 2 . From the data, empirical exponential forms 

V »=i 

r,T wlO Q+/3h , aw -0.40, /3 « 0.47 (6) 

can be approximately formulated [the solid curve in Fig.l]. It is emphasized that 
the period T"s (also the transient time r's) are much shorter than those allowed by 
the spaces of the computer data for the given precisions. For instance, we have T 
« 2 x 10 3 << 10 8 for 10~ 8 precision and T w 10 6 << 10 13 for 10~ 13 precision. 




5 10 15 

h 



Fig. 1. Average transient time r's (circle) and average period T's (square) for entering periodicity 
of numerical realizations of chaotic trajectories of Eq.(3) for different precisions. Solid line is the 
result of the empirical formulas (6), which fits the numerical data satisfactorily. Vertical bars are 
errors of statistical data of average r and T for different precisions. 

There are several interesting observations in Fig.l and Eq.(6) worthwhile re- 
marking. First, the average period length has a power law scaling against the dis- 
crete size e = 10 - ' 1 , as T oc e _/3 , (3 w —-^(within the statistical error ±0.02) 
with Z?2 being the correlation dimension of the given chaotic attractor (in our case 
D 2 = 1), this agrees with the theoretical prediction of [8,11]. Moreover, the average 
transient time r has the similar scaling as T. From the point of view of cryptography, 
transient time r has the same significance as period T itself. For a given discrete 
size 10 , there arc a number of periodic orbits of which the corresponding periods 
vary largely. There exist some orbits with very small periods even h is relatively 
large. Large r is useful for keeping the communication security even in the case of 
small T, since it is not the quantity T but the quantity T + r that is related to the 
security problem of the chaotic cryptography. 
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Another interesting feature of CRCT is that the number of periodic orbits which 
the system eventually enters for a given computation precision is very small. An 
intuitive picture is: Since there are infinitely many unstable periodic orbits densely 
embedded in any chaotic attractor it is easy for a trajectory to be trapped in one 
of these orbits after a course grain caused by finite precision computation, and thus 
the number of realized periodic orbits may be huge as the computation precision 
is high (h 3> 1). To our great surprise, for any precisions (up to the precision 
10~ 14 ) the numbers of observed distinctive orbits are very small, and for most of 
truncation precisions we observe that one orbit of all periods takes major part of 
probability, and dominates the statistical behavior of CRCT. These results are listed 
in Table 1, which are strangely against the above intuition. Note, the existence of 
such a dominant period dose not support the theoretical results on the probability 
of periods in [8,10]. The reason for this disagreement may be the following. In [8] 
the prediction is proven based on approximating the map to be deterministic in 
time and purely random in variable space map (though deterministic [8,10]). The 
existence of dominant period breaks the validity of this intuitive picture. 



Table 1. Period lengths and corresponding probabilities by different decimal precisions 
realization. 



Precision 


Number 


Period length Ti (ratio) 


e = 


w- h 


of tests 




10" 


-1 


10 1 


1(1.00) 


10 


-2 


10 2 


3(0.60) 1(0.40) 


io- 


-3 


10 3 


1(0.92) 13(0.08) 


io- 


-4 


10 4 


103(0.4916) 97(0.3298) 1(0.1760) 6(0.0014) 2(0.0012) 


io- 


-5 


10 5 


227(0.92279) 1(0.06093) 29(0.00937) 
23(0.00661) 5(0.0016) 4(0.00008) 2(0.00006) 


io- 


-6 


10 6 


155(0.680588) 1(0.210002) 16(0.065585) 88(0.033479) 
79(0.00826) 8(0.001122) 30(0.000512) 11(0.000336) 6(0.000116) 


10" 


-7 


10 3 


1078(0.505) 136(0.427) 1(0.068) 


io- 


-8 


10 3 


2412(0.756) 419(0.122) 1(0.10) 428(0.018) 118(0.004) 


10 


-9 


10 3 


5957(0.97) 1(0.028) 860(0.002) 


10~ 


10 


500 


26358(0.814) 1319(0.132) 1680(0.026) 4317(0.012) 
1(0.012) 1233(0.002) 643(0.002) 


10~ 


11 


200 


101487(0.96) 22916(0.035) 2982(0.005) 


io- 


12 


200 


1(0.815) 95914(0.185) 


io- 


13 


66 


1(0.818) 985982(0.182) 


io- 


14 


28 


2385320(0.429) 1908178(0.357) 183458(0.107) 1(0.107) 



In the above discussion we used the decimal representation and the fixed- 
precision truncations. The actual computer round-off discretization use the binary 
representation and relative-precision truncations. In Fig. 2 we do the same as Fig.l 
by applying binary representation and relative-precision truncations. The results are 
similar to those of Fig. 1 , showing that though the discretization of variable space 
has some significant changes in the dynamics and statistics of chaotic systems, the 
detailed methods of discretization is not crucial. 



3, 2008 16:25 WSPC/INSTRUCTION FILE ctextemparxiv 



6 S.H. Wang 




Fig. 2. The same as Fig.l with binary representation and relative-precision truncations applied. 
The statistical results are unchanged from Fig.l. 

3. Periodicity of single chaotic maps in double precision computer 
realization 

Wo directly and numerically compute Eq.(l) in Table 2 by the double precision 
2~ 52 computations, and by using the arithmetic of [20] to detect the period length. 
This method cannot exactly give transient time, but can record period exactly 
and effectively. A significant advantage of this approach is that it can compute 
extremely long period with small storage space and fast computation. We can see 
only 12 different periods for 10 5 different initial conditions from Table 2. There 
exists a major probability period and there are near to zero probabilities for some 
periods. The average period length is about 5.7 x 10 6 iterations, which fits form (6) 
approximately. The results of the direct double precision computing well confirm 
our discretization method. We use C code to validate it. It is interesting to point 
out that for the true double precision computation of the model Eq.(l), the period 
of length 5638349 (which is one order smaller than the period T ~ 2 26 « 6 x 10 7 
predicted by [8]) can appear with probability 67.75% (10 5 different initial conditions 
are computed). This knowledge may be useful for scientists interested in computer 
simulations of chaotic systems. We also consider the single tent map s n+ i = 1 — 
|1.95a;„ — 1| as our example and compute the period length in Table 2 by the double 
precision. There are the similar results as the single logistic map. 

4. Periodicity of chaotic coupled map lattices in Computer 
realization 

In Figs.l and 2 it is clearly observed that with a single chaotic map the period of 
orbits of computer realization is rather short. For double precision (2~ 52 ) the average 
period is about of order 10 6 ~ 7 , which can be very easily reached. This short period 
problem may cause serious disadvantages in chaos applications. For instance, a single 
autonomous chaotic map can never allow high security in computer realization 
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Table 2. Period lengths and corresponding numbers of 10 5 different initial 
conditions of logistic and tent maps by double precision realization. 



T 


IN umbers 


IN umbers 


NumbersNumbers 


IN Ullll KTS 


(Logistic rricip) 


nf 9nnnn 


onnnn 

Ol zuuuu 


of 20000 of 20000 


nf 9nnno 


5638349 


13533 


13672 


13533 


13482 


13631 


1 


3233 


3312 


OZoO 


oolo 


3253 


14632801 


2339 


2297 


2284 


2272 


2172 


10210156 


336 


346 


349 


346 


357 


2441806 


284 


290 


302 


285 


275 


2625633 


264 


274 


285 


290 


301 


420909 


7 


1 


3 


3 


3 


960057 


4 


3 


5 


5 


6 


1311627 





4 


2 


2 


2 


4389 





1 


1 








510250 








1 


1 





234209 











1 





Ti 


Numbers 


Numbers 


Numbers 


Numbers 


Numbers 


(Tent map) 


of 20000 


of 20000 


of 20000 


of 20000 


of 20000 


216135802 


10111 


10105 


10089 


10045 


10053 


1009008 


7965 


7966 


7968 


8006 


7993 


128662357 


1820 


1824 


1836 


1847 


1851 


8703769 


47 


48 


49 


45 


44 


15083761 


35 


34 


36 


35 


38 


4387781 


15 


16 


15 


15 


15 


3404560 


6 


6 


6 


6 


6 


1299622 


1 


1 


1 


1 






of chaotic cryptography since long period is a necessary condition for any kinds 
of cryptosystems [19-21]. In Ref.22, the authors suggested to use spatiotemporal 
chaotic systems for cryptography and declared that they can reach very high level 
of security, it is then interesting to investigate the periodicity of CRCTs of such 
spatiotemporal systems. 

Let us extend the single map (1) to one-way coupled map lattice 

x n+1 (i) = (l-e)f[x n (i)]+ef[x n (i-l)}, i = l,2,...,JV (7) 
f(x) = 4x(l-x), x(0)=x(N) 

where we use periodic boundary condition of system size N , and take e = 0.95. 
This system is used for chaotic encryption in [22]. 

We are now interested in how both the system size N and the computer precision 
10- h affect the period T and the transient time r of CRCTs. In Figs.3(a)-(c) we 
do the same as Fig.l with N = 2, N = 3 and N — 4 in Eq.(7), respectively. The 
average lengths of period and the transient iterations increase exponentially with 
h much faster than Fig.l. Moreover, for the same precision, the average lengths of 
the period and the transient are longer if the system size is taken larger. In Fig. 3(d) 
we fix h = 3 and plot T and r vs N, an exponential increasing of T and r with N 
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are confirmed. In all figures of Fig.l and Fig. 3, the solid curves are drawn from the 
united empirical form 

rW.Tjf 1 oc lO^+W" = 10 QAr e-^, aw -0.40, « 0.47 (8) 

which coincides satisfactorily with the numerical data for wide range of N and h. 
By increasing the system size N, the coupled systems can manifest its chaoticity in 
much longer time though periodicity must appear for sufficiently large time scale. If 
the system size is sufficiently large the periodic behavior of CRCTs is practically not 
observable. According to the empirical formula (8), for N = 4, and double precision 
e = 2~ 52 , the periodicity occurs after 10 28 iterations, and a common PC (with 1.7G 
CPU) needs 10 12 years, to make so many iterations. Thus, the short period problem 
of CRCTs can be practically solved by applying coupled (spatiotemporal) chaotic 
systems with a sufficiently large N. This is one of several most important reasons 
why the system in [22] may reach very high security. 




Fig. 3. (a)-(c) The same as Fig.l with Eq.(7) considered, s = 0.95, AT = 2, N = 3, and N = 4, 
respectively, (d) T and r plotted vs N for Eq.(7) with h = 3. Solid lines are the results of the 
empirical formulas (8), which fit well the numerical data. It is shown that the average transient r 
following the power law Eq. (8) satisfactorily with smaller fluctuation than that of T. 

The phenomenon of small number of periodic orbits of CRCT is also observed 
for coupled maps. In Fig. 3 for N > 1 and relatively large h, we can try only small 



3, 2008 16:25 WSPC/INSTRUCTION FILE ctextemparxiv 



Instructions for Typing Manuscripts (Paper's Title) 9 

number of tests for each given precision h. Therefore, the fluctuation of T is rather 
larger. Sometimes we obtain very small T's, which deviate from the exponential law 
(the solid line of Fig. 3) considerably. However, almost in all our tests the transient 
time t's follow Eq.(8) satisfactorily with fluctuation much smaller than that of T's 
(see data of Fig. 3). Since the cryptography of chaotic systems is related to T + r 
rather than T, such nice power law behavior of r is useful for controlling the security 
of cryptosystems. 

5. discussion and conclusion 

In conclusion we have investigated the problems of periodicity and transient to pe- 
riodicity of chaotic trajectories in computer realization with finite computation pre- 
cision. It is found that both average period and the transient time to periodic orbits 
have power law relation with the computer round-off precision. For low-dimensional 
systems the periods of chaotic trajectories and the corresponding transients can be 
rather short even the double-precision computation is applied. And for each pre- 
cision the number of periodic orbits of computer realization is small even if the 
precision is rather high and the number of tests is very large. The problems of short 
period and small number of periodic orbits may seriously affect the applications of 
chaos, e.g., it results in weak cryptography in secure communication. Nevertheless, 
this problem can be satisfactorily solved by applying spatiotemporal chaos (or say, 
coupled chaotic systems). With few coupled chaotic subunits the period of CRCTs 
may be as large as practically unreachable. A surprising observation worthwhile 
remarking is that by applying different initial conditions there always exists a dom- 
inant period (for all tested computation precisions and sizes of coupled systems) 
which appears with major portion of probability. 

This work was supported by the National Science Foundation of China and 
Nonlinear Science Project. 

References 

1. A.N. Kolmogorov, D. Akad and Nauk, USSR 98, 1 (1954). 

2. E.N. Lorenz, J. Atmos. Sci. 20, 130 (1963). 

3. D. Ruelle and F. Takens, Commun. Math. Phys. 20, 167 (1971). 

4. M.J. Feigenbaum, J. Stat. Phys. 19, 25 (1978). 

5. J.P Eckmann and D. Ruelle, Rev. Mod. Phys. 57, 617 (1985). 

6. J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems and Bi- 
furcations of Vector Fields, Springer- Verlay, (1983). 

7. B.L. Hao and W.M. Zheng, Applied Symbolic Dynamics and Chaos, World Scientific, 
Singapore (1998). 

8. C. Grebogi, E. Ott, and A. Yorke, Phys. Rev. A 38, 3688 (1988). 

9. C. Beck, Phys. Lett. A 136, 121 (1989). 

10. U. Bastolla, and G. Parisi, J. Phys. A: Math. Gen., 30 (1997). 

11. C. Dellago, and W.G. Hoover, Phys. Rev. E 62, 6275 (2000). 

12. L.M. Pecora and T.L. Carroll, Phys. Rev. Lett. 64, 821 (1990). 

13. L.M. Cuomo and A.V. Oppenheim, Phys. Rev. Lett. 71, 65 (1993). 



3, 2008 16:25 WSPC/INSTRUCTION FILE ctextemparxiv 



10 S.H. Wang 

14. L. Kocarev, K.S. Halle, K. Eckert, L.O. Chua and U. Parlitz, Int. J. Bif. and chaos 
2(3), 709 (1992). 

15. L. Kocarev and U. Parlitz, Phys. Rev. Lett. 74(25), 5028 (1995). 

16. J.H. Xiao, G. Hu and Zh.L. Qu, Phys. Rev. Lett. 77, 4162 (1996). 

17. D.G. Van Wiggeren and R. Roy, Science 279(20), 1198 (1998). 

18. S. Sundar and A. A. Minai, Phys. Rev. Lett. 85(25), 5456 (2000). 

19. R. Matthews, Cryptologia 13, 29 (1989) . 

20. D. Wheeler, Cryptologia 13, 243 (1989) . 

21. F. Dachselt, W. Schwarz, IEEE trans. Circuits syst. I, 48(12), 1498 (2001). 

22. S.H. Wang, J.Y. Kuang, J.H. Li, Y.L. Luo, HP. Lu and G. Hu, Phys. Rev. E 66, 
065202 (2002) 



